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ABSTRACT 

A method for extracting the Levy stability index /i from the multi-fractal 
spectrum f(a) in high energy multiparticle production is proposed. This index 
is an important parameter, characterizing the non-linear behaviour of dynamical 
fluctuations in high energy collisions. Using the random cascading a model as 
example, the validity of this method is tested. It is shown that this method, 
basing on a linear fit, is consistent with and more accurate than the usual method 
of fitting the ratio of qth to 2nd order multi-fractal (Renyi) dimensions to the 
Peschanski formula. 
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Recently, the search, starting from eighties' 1 !, f or the non-linear phenomena — inter- 
mittency or fractality in the multiparticle final states of high energy collisions' 2 ' has had 
breakthrough' 3 '. It is found that dynamical fluctuations do exist in high energy hadron- 
hadron collisions, but is aniosotropic rather than isotropic' 4 !. Power law in the higher- 
dimensional factorial moment versus partition number has been observed in hadron-hadron 
collisions, when the anisotropy of dynamical fluctuations is taken into account properly' 5 ' 6 !. 
On the other hand, the dynamical fluctuations in e + e _ collisions are found to be approxi- 
mately isotropic ' 7 1 . 

As a (self-affine or self-similar) fractal system, the multiparticle final state in high en- 
ergy collisions can be characterised by an important parameter — the Levy stability index 
/i' 8 l . This parameter tells us the behaviour of elementary fluctuations at the tail of distri- 
bution. To extract its value from exprimcntal data more reliably is an essential task for the 
understanding of fluctuation dynamics. 

Previously, people usually use the anomalous scaling of factorial moment' 9 ' 

j? _ 1 (n m {n m - 1) ■ ■ ■ (n m - q + 1)) 

M 2-/ {nm)q W 

to get the scaling index (IM index) (j) q and multi-fractal (Renyi) dimensions^ 10 ! 

D " = 1 -^h (2) 

of integer order q, and then fit the ratio (1 — D q )/(1 — D 2 ) to the formula' 11 ' 



1-D q = 1 g"-g 
1 - D 2 q - 1 2** - 2 



(3) 



to get the value of /i. In practice only a few indices <f> q of integer q > 2 (q = 3, 4, 5 e.g.) can 
be obtained and the resulting value of /i basing on a fit of very few points to the complicated 
formula (3) is obviously unsatisfactory. 

In order to establish an alternative way for extracting the value of fi, let us notice that 
the gth order multi-fractal Renyi dimension D q of any order q can be obtained from the 
multifractal spectrum /(a) through' 10 ' 

D q = -^j(qa-f(a)). (4) 

Therefore, it should be able to get the Levy stability index /x directly from the multifractal 
spectrum /(a). This short note is aimed to discuss this problem. 
It can be seen from Eq.(4) that 

r q = (q- l)D q (5) 



is the Legandre transform of f(a). Therefore, we have' 12 ' 



f(a) = qa q - T q . (7) 

From Eq.(3) we have 

, — q 
q - 1 

where A = (1 — D 2 )/(2 fJ - — 2). Inserting into Eq.(5) and making use of Eq's. (6), (7) we get 

r q = -! + {! + A)q-Aq^ (8) 

a q = 1 + A-iiAq*- 1 , (9) 
/(a) = l-( M -W. (10) 



Therefore, defining 
we finally get 



B=l+A=l+- — — , (11) 
2^-2 v ; 



1 - f(a) oc (B - a )f/(f-i) when a < B. (12) 



Thus we propose the following procedure for getting the Levy stability index \x from the 
multi- fractal spectrum f(a) : 

(1) Find out the value of a where /(a) = 1, call it B . 

(2) Fit the a < B part of ln(l — f(a)) versus ln(_B — a) to a straight line and get the 
slope C. 

(3) Get Levy stability index as 

As an example, we apply this method to the random cascading a-modell 13 ! . This model 
describes each multiparticle event as a series of steps, in which the initial phase space region 
A is repeatedly divided into A = 2 parts. After v steps we get M — 2 V sub-cells of size 
5 = A/M. At each step s the normalized particle density is obtained in each of the two 
parts by multiplication of the normalized density in the step s — 1 by a particular value of 
the random variable w„j„, where j v is the position of a sub-cell at the vth step (1< j v < 
2 U ). The elementary fluctuation probability uj is chosen as' 13 ' 

w„, 2 j-i - ^(1 + ar) ; Wv,2j = \0-~ ar )> ( 14 ) 

in which, r is a random number distributed uniformly in the interval [—1, 1], a is a model 
parameter, characterizing the strength of fluctuations. In our calculation, we have chosen 
its value as a = 0.3. 

The resulting multi-fractal spectrum /(a) in this model are shown in Fig. 1. The value 
of a for f(a) = 1 is found to be B = 1.0113. The fitting of (1 - D q )/(1 - D 2 ) to Eq.(3) and 
the linear fit of ln(l — f(a)) versus ln(£? — a) are shown in Fig's. 2 and 3 respectively. The 
resulting Levy indices [i from these two methods with different regions of q together with 
the corresponding % 2 /DF are 



Q 



Fit from (1 - D q )/(1 - D 2 ) 

{q = 2,3,4,5) fi = 1.929 ± 0.059 x 2 /DF = 0.016 /2. 

Linear fit of ln(l — f(a)) versus ln(_B — a) 

(q = 0.1-5, a = 1.009-0.909) (j, = 1.962 ± 0.013 x 2 /DF = 0.758 /48. 

(q = 2-5, a = 0.967-0.909) /z = 1.950 ± 0.022 xV DF = 0.107 /29. 

(q = 3-5, a = 0.945-0.909) (j, = 1.940 ± 0.015 x 2 /DF = 0.015 /19. 
cf. Fig. 4. 

Thus we see that the method we proposed for extracting the Levy stability index [i from 
the multi- fractal spectrum f(a), basing on a linear fit, is consistent with the usual method 
of fitting the ratio of gth to 2nd order multi-fractal (Rcnyi) dimensions to the Peschanski 
formula Eq.(3). It is more accurate and reliable and is reasonably stable with respect to the 
region of q or a used in fitting. 

This method can readily be applied to real experimental data analysis as a better way 
for getting Levy stability index. 
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Figure Captions 

Fig.l Multi- fractal spectrum 

Fig.2 Fitting of (1 -D q )/ (1 - D 2 ) to Eq.(3) for q = 2, 3, 4, 5 
Fig.3 Linear fit of ln(l - f(a)) vs. ln(B - a) for q = 0.1 - 5 

Fig. 4 The resulting \x value from different ways of fitting together with the q region used 
and the x 2 / DF 




Fig.3 Linear fit of ln(l - /(a)) vs. ln(B - a) for q = 0.1 - 5 



Fit from Region of q value 
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Fig. 4 The resulting fj, value from different ways of fitting 
together with the q region used and the x 2 /DF 
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